Thermodynamic Properties of Kagome Lattice in ZnCu3(OH) 6 Cl2 Herbertsmithite 
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Strongly correlated Fermi systems are among the most intriguing and fundamental systems in 
physics, whose realization in some compounds is still to be discovered. We show that herbertsmithite 
ZnCu3(OH)6Cl2 can be viewed as a strongly correlated Fermi system whose low temperature ther- 
modynamic in magnetic fields is defined by a quantum critical spin liquid. Our calculations of its 
thermodynamic properties are in good agreement with recent experimental facts and allow us to 
reveal their scaling behavior which strongly resembles that observed in HF metals and 2D 3 He. 
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An explanation of the rich behavior of strongly cor- 
related Fermi systems still continues to be among the 
main problems of the condensed matter physics. One of 
the most interesting and puzzling issues in the research of 
strongly correlated Fermi systems is the non- Fermi liquid 
(NFL) behavior detected in their thermodynamic prop- 
erties. Under the application of external fields, e.g. mag- 
netic field B, the system can be driven to a Landau Fermi 
liquid behavior (LFL). Such a behavior was observed in 
quite different objects such as heavy-fermion (HF) metals 
and two-dimensional 3 He [2H5fl 

Recently the herbertsmithite ZnCu3(OH)eCl2 has 
been exposed as a S = 1/2 kagome antiferromagnet @ 
and new experimental investigations have revealed its un- 
usual behavior Because of the electrostatic envi- 
ronment, Cu 2+ is expected to occupy the distorted oc- 
tahedral kagome sites. Magnetic kagome planes Cu 2+ 
5 = 1/2 are separated by nonmagnetic Zn 2+ layers. Ob- 
servations have found no evidence of long range magnetic 
order or spin freezing down to temperature of 50 mK, in- 
dicating that ZnCu3(OH)fiCl2 is the best model found of 
quantum kagome lattice [7H10j . The specific heat C, aris- 
ing from the Cu spin system, at T < 1 K appears to be 
governed by a power law with an exponent which is less 
than or equal to 1. At the lowest explored temperature, 
namely over the temperature range 106 < T < 400 mK, 
C follows a linear law temperature dependence, C oc T, 
and for temperatures of a few Kelvin and higher, the spe- 
cific heat becomes C(T) oc T 3 and is dominated by the 
lattice contribution [7H9(. At low temperatures T < 1, 
the strong magnetic field dependence of the specific heat 
C suggests that C is predominately magnetic in origin 
There are a number of papers suggesting that the 
5=1/2 model on the kagome lattice can be viewed 
as the gapless critical spin liquid 0-i, EMI- These 
facts allow us to test both the NFL and LFL behavior of 
ZnCu3(OH)gCl2 and to show that a Fermi quantum spin 
liquid formed in the herbertsmithite determines its low 
temperature thermodynamic properties. 

Contrary to the C oc T behavior 0, S|, the ob- 
served spin liquids contribute a T 2 specific heat which 



in the model is not sensitive to an applied magnetic field 
[l5L llil ]. Moreover the magnetic susceptibility x(T) of 
ZnCu3(OH)gCl2 shown in Fig. [1] displays an unusual be- 
havior [9(. At B > 3 T, x(T) has a maximum Xmax(r) 
at some temperature T max (i3). The maximum Xmax(^) 
decreases as magnetic field B grows, while T max (i3) shifts 
to higher T reaching 15 K at B = 14 T. At B < 1 T as 
seen from Fig. HJ y{T)<x T~ a with a = 2/3. The cal- 
culated exponent [2|, Il7| is in good agreement with the 
experimental value a = 2/3 ~ 0.66 [9(. The observed 
behavior of x strongly resembles that in HF metals and 
is associated with their proximity to a quantum critical 
point (QCP) 0, IT7I fljj |. As a result, we safely assume 
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FIG. 1: Temperature dependence of the magnetic suscepti- 
bility x a t different magnetic fields for ZnCu3(OH)gCl2 0. 
The illustrative values of Xmax and Tmax at B = 7 T are also 
shown. Our calculations made at B — are depicted by the 
solid curve representing x(T) oc T~ a with a — 2/3. 

that a quantum deconfined Fermi critical spin liquid with 
essentially gapless excitations formed by neutral fermions 
is realized in ZnCu3(OH)6Ci2 and located very near QCP 
0. Thus, ZnCu 3 (OH) 6CI2 turns out to be located at its 
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QCP without tuning this substance to QCP using a con- 
trol parameter such as magnetic held, pressure, or chem- 
ical composition. This observation is in sharp contrast to 
the common practice applied to tune HF metals to their 
QCP's. A simple kagome lattice may have a dispersion- 
less topologically protected branch of the spectrum with 
zero excitation energy that is the flat band jll, 12 1. In 
that case a fermion condensation quantum phase tran- 
sition (FCQPT) 0] can be considered as QCP of the 
ZnCu3(OH)gCl2 quantum spin liquid. 

In this letter we uncover the quantum spin liquid phase 
and its QCP in the herbertsmithite ZnCu3(OH)gCi2 and 
explain its low temperature thermodynamic in magnetic 
fields. For the first time we calculate the susceptibility 
X, magnetization M and specific heat C as functions of 
temperature T versus magnetic field B. Our calculations 
are in good agreement with the experimental facts and 
allow us to reveal their scaling behavior which strongly 
resembles that observed in HF metals and 2D 3 Hc. 

To study the low temperature thermodynamic and 
scaling behavior, we use the model of homogeneous 
heavy- fermion liquid Q. This model permits to avoid 
complications associated with the crystalline anisotropy 
of solids. We propose that the critical quantum spin liq- 
uid is composed of fermions. These fcrmions with zero 
charge and spin a = 1/2 occupy the corresponding Fermi 
sphere with the Fermi momentum pp. The ground state 
energy E(n) is given by the Landau functional depending 
on the quasiparticle distribution function n a (p), where p 
is the momentum. Near FCQPT point, the effective mass 
M* is governed by the Landau equation 0, Ell 
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Here we have rewritten the quasiparticle distribution 
function as n a (p, T, B) = n a (p,T = 0, B = 0) + 
5n a (p, T, B). The Landau amplitude F is completely 
defined by the fact that the system has to be at QCP 
of FCQPT @, 03, ED]. The sole rolc of the Landau am- 
plitude is to bring the system to FCQPT point, where 
Fermi surface alters its topology so that the effective mass 
acquires temperature and field dependence [H, E3, [HI- 
At this point, the term 1/M* vanishes and Eq. (JlJ be- 
comes homogeneous. It can then be solved analytically 
@,[I3l- At _B = 0, the effective mass str ong ly depends on 
T demonstrating the NFL behavior 0, 1 1 7t | 



M*(T) ~ a T T 



-2/3 
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At finite T, the application of magnetic field B drives the 
system the to LFL region with 



At finite B and T near FCQPT, the solutions of Eq. 
([T]) M*{B,T) can be well approximated by a simple uni- 
versal interpolating function. The interpolation occurs 
between the LFL (M*(T) oc const) and NFL (M*(T) oc 
T -2 / 3 ) regions @, E3]- ^ i s convenient to introduce the 
normalized effective mass M^ and the normalized tem- 
perature T/v dividing the effective mass M* by its max- 
imal values, M M , and temperature T by T max at which 
the maximum occurs. Equation ([T]) allows us to calculate 
the thermodynamic properties for the normalized suscep- 



Sincc C/T oc M* , the 
We note that our cal- 



tibility xn = x/Xmax = Ab- 
normalized (C/T) = xn = M^ 
culations of M^ based on Eq. ([I]) do not contain any 
free fitting parameters. The normalized effective mass 
Mjy = M*/M M as a function of the normalized temper- 



ature y = T N = T/T max is given by @, [n| 



M* N (y) « c - 



ciy 



c-iy 



8/3 ' 



(4) 



Here Cq = (1 + c%)/ (1 + C\), C\ and c-i are fitting param- 
eters, approximating the Landau amplitude. Magnetic 
field B enters Eq. ([T]) only in the combination ubB/IcbT, 
making fcsT max ~ fisB where ks is the Boltzmann con- 
stant and fib is the Bohr magneton @, [T3]- Thus, in the 
presence of magnetic fields the variable y becomes 



y = t/t u 



k B T/fi B B. 



(5) 



M*(B) ~ a B B 
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The variables T and B enter Eq. (J5J) symmetrically; 
therefore Eq. @ is valid for y = fi B B /k B T . In what 
follows we use Eq. ((4]) to clarify our calculations based 
on Eq. |(1]). It follows directly from Eqs. (gj and 

© that x(k B T/fi B B)T 2 / 3 oc y 2 ' z M* N (y). Since the 
magnetization M(B,T) = J x{B,T)dB, we obtain that 
M(B,T)T -1 / 3 depends on the only variable y. These 
observations confirm the scaling behavior of both xT 66 
and MT -0 - 34 that is experimentally established in Q. 

We are now in a position to construct the schematic 
phase diagram of ZnCu3(OH)gCl2. The phase diagram 
is reported in Fig. At T = and B = the sys- 
tem is located at QCP of FCQPT without tuning. Both 
magnetic field B and temperature T play the role of the 
control parameters, shifting the system from its QCP and 
driving it from the NFL to LFL regions as shown by the 
vertical and horizontal arrows. At fixed temperatures the 
increase of B drives the system along the horizontal ar- 
row from the NFL region to LFL one. On the contrary, 
at fixed magnetic field and increasing temperatures the 
system transits along the vertical arrow from the LFL 
region to NFL one. The inset to Fig. [2] demonstrates the 
behavior of normalized effective mass M^j versus normal- 
ized temperature Tjy following from Eq. ((4]). It is seen 
that the temperature region T/v ~ 1 represents a transi- 
tion region between the LFL behavior with almost con- 
stant effective mass and the NFL one, having the T~ 2 / 3 
dependence. It is seen from Eqs. (0| and |5]) and Fig. [5] 
that the width of the transition region T w oc T oc B. 
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ior. Consistent with the phase diagram displayed in Fig. 
[2] at growing temperatures (y ~ 1) the LFL behavior 
first converts into the transition one and then disrupts 
into the NFL behavior. This demonstrates that the spin 
liquid of ZnCu 3 (OH) 6 Cl2 is close to QCP and behaves 
as the HF liquid of YbRli2Si2. It is seen, that the low- 
temperature ends (T/v < 0.5) of the curves do not merge 
and their values decrease as B grows representing the full 
spin polarization of the HF band at the highest reached 
magnetic fields (H)]. Indeed, at low Tjy, xn at B = 14 T 
is close to (C/T)n at B = 18 T, while our calculations 
shown by the solid curve is close to both functions. 



FIG. 2: Schematic phase diagram of ZnCu3(OH)eCl2. The 
vertical and horizontal arrows show LFL-NFL and NFL-LFL 
transitions at fixed B and T respectively. Inset shows a 
schematic plot of the normalized effective mass versus the nor- 
malized temperature. Transition region, where reaches 
its maximum at Tjv = T/T max = 1, is shown by the arrows 
and hatched area in both the main panel and in the inset. 
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FIG. 4: Normalized magnetization Mjv(y) collected on mea- 
surements on ZnCu 3 (OH) 6 Cl2 0| and YbRh 2 Si 2 [H] at differ- 
ent temperatures shown in the corresponding legends. Shown 
by the arrow a kink is clearly seen at y ~ 1. The normalized 
entropy S'jv(y) is extracted from measurements on 2D 3 He @] 
at different densities x shown in the legend. The solid curve 
represents our calculations of the normalized magnetization. 



FIG. 3: Normalized susceptibility \n = x/Xmax = M% 
versus normalized temperature Tjv. xn is extracted from 
the measurements of the magnetic susceptibility x i n mag- 
netic fields B [sj shown in Fig. [TJ Normalized specific heat 
(C/T)jv = M^f is extracted from the measurements of C/T on 
YbRh2Si2 in magnetic fields B [2^. The corresponding fields 
B are listed in the legends. Our calculations made at field B 
completely polarizing the quasiparticle band are depicted by 
the solid curve tracing the scaling behavior of M^. 

The experimental data on measurements of xn @, 
(C/T) w = M$j [20t and our calculations of M* N at fixed 
magnetic field B that completely polarizes the quasipar- 
ticle band are shown respectively by the geometrical fig- 
ures and solid curve in Fig. [3] It is clearly seen that the 
data collected on both ZnCu3(OH)6Cl2 and YbRli2Si2 
collapse into the same curve, obeying the scaling behav- 



Both the normalized magnetization Mjy(y) = 
M(B/Bk)/M(Bk), extracted from measurements of the 
magnetization M(B) [§] depicted by the geometrical fig- 
ures, and calculated M^{y) shown by the solid line are 
reported in Fig. g] Here, Tk is the temperature at which 
the magnetization demonstrates the kink, while the sys- 
tem enters the transition region shown in Fig. CU [2]. The 
normalized entropy Sn(u) — S(T/Ti n f)/S(Ti n f) is ob- 
tained from measurements of the entropy S on 2D 3 He 
Here T in j is the temperature at which the system 
enters the transition region and S possesses its inflection 
point, clearly seen in the data (see the Supporting Online 
Material for Fig. S8, A). It is seen from Eq. (gj and 
from the inset to Fig. [2l that at y < 1, Sn = Mjv oc y, 
and at y > 1, Sn = Mn oc y 1 ^ 3 - This behavior produces 
the kink and makes the scaled data merge into a single 
curve in terms of the variable y Q ■ Our calculations are 
in good agreement with the measurements. 
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FIG. 5: Panel A: The temperatures T max (_B) at which the 
maxima of \ ( see Fig- are located. The solid line represents 
the function r max oc aB, a is a fitting parameter, see Eq. 
([5]). Panel B: The maxima Xmax of the function x(T) versus 
magnetic field B (see Fig. [TJ). The solid curve is approximated 
by Xmax(S) = dB~ 2 / 3 , see Eq. (J3|, d is a fitting parameter. 
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FIG. 6: The specific heat C(B, T) versus magnetic field B 
measured on ZnCu3(OH)eCl2 at two different temperatures T 
0] listed in the legends is shown by the triangles and squares. 
Our calculations are depicted by the solid curves tracing the 
LFL behavior of C(B,T) = ai_B" 2/3 T, see Eq. ©, with oi 
being the fitting parameter. 



In Fig. [5j panel A, the solid squares denote tempera- 
tures T , max (B) at which the maxima of x(^) occur ver- 
sus magnetic field B. In the panel B, the corresponding 
values of the maxima Xmax(-B) are shown by the solid 
diamonds versus B. It is seen that the agreement be- 
tween the theory and experiment is good in the entire 
magnetic field domain. Our calculations of the specific 
heat C(B, T) are shown in Fig. [BJ For T of a few Kelvin 
and higher, the lattice contribution to the specific heat 



is the most significant contribution. However, this con- 
tribution diminishes at low temperatures, and at T < 1 
K C is predominately formed by the spin liquid Q . It is 
seen from Fig. [1JJ that in the LFL region at fc^T < figB, 
C{B,T) oc M*T oc B- 2 / 3 T, and field B completely de- 
fines the M*{B) behavior given by Eq. ([3]). Clearly, our 
calculations are in good agreement with the measure- 
ments when the system demonstrates the LFL behavior. 
Indeed, at T = I K the system exhibits the LFL behav- 
ior at B > 2 T, while at T = 0.5 K the LFL behavior is 
observed even at lower values of B, namely B > 1 T. 

In summary, we have shown that the kagome lattice 
of ZnCu3(OH)@Cl2 can be viewed as a strongly cor- 
related Fermi system whose thermodynamic is defined 
by the quantum critical spin liquid located at FCQPT. 
Our calculations of the thermodynamic properties are in 
good agreement with the experimental facts and their 
scaling behavior coincides with that observed in HF 
metals and 2D 3 He. We have also demonstrated that 
ZnCu 3 (OH) 6 Cl 2 exhibits the LFL, NFL and the transi- 
tion behavior as HF metals and 2D 3 Hc do. 
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